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ABSTRACT 



Context. As a galaxy evolves, its chemical composition changes and the abundance ratios of different elements are powerful probes 

of the underlying evolutionary processes. Phosphorous is an element whose evolution has remained quite elusive until now, because 

it is difficult to detect in cool stars. The infrared weak Pi lines of the multiplet 1, at 1050-1082 nm, are the most reliable indicators of 

the presence of phosphorus. The availability of CRIRES at VLT has permitted access to this wavelength range in stellar spectra. 

Aims. We attempt to measure the phosphorus abundance of twenty cool stars in the Galactic disk. 

Methods. The spectra are analysed with one-dimensional model-atmospheres computed in Local Thermodynamic Equilibrium (LTE). 

The line formation computations are performed assuming LTE. 

Results. The ratio of phosphorus to iron behaves similarly to sulphur, increasing towards lower metallicity stars. Its ratio with respect 

to sulphur is roughly constant and slightly larger than solar, [P/S]=0. 10 ± 0.10. 

Conclusions. We succeed in taking an important step towards the understanding of the chemical evolution of phosphorus in the 

Galaxy. However, the observed rise in the P/Fe abundance ratio is steeper than predicted by Galactic chemical evolution model model 

developed by Kobayashi and collaborators. Phosphorus appears to evolve differently from the light odd-Z elements sodium and 

aluminium. The constant value of [P/S] with metallicity implies that P production is insensitive to the neutron excess, thus processes 

other than neutron captures operate. We suggest that proton captures on ^^"Si and a captures on ^^Al are possibilities to investigate. 

We see no clear distinction between our results for stars with planets and stars without any detected planet. 

Key words. Stars: abundances - Stars: atmospheres - Line: formation - Galaxy: evolution - Galaxy: disk - radiative transfer 



1. Introduction 

Phosphorus is an abundant element in the Universe, and in the 
solar photosphere it is among the top twenty most abundant ele- 
ments. In the periodic table of elements, phosphorus, with Z=15, 
is in the same group as nitrogen, lying between silicon (Z=14) 
and sulphur (Z-16). Silicon is a well-studied element, and its 
abundance is easily measured in photospheres of stars of type 
F, G, or K, by analysing the Sii lines. In the last fifteen years, 
sulphur has been systematically investigated by several groups 
that analysed the few multiplet s of S i available in the observed 
spectra (for further details see ISpite et al]|2011l and references 
therein). This is not the case for phosphorus, that, before this 
work, had never been analysed sys tematically in cool stars. The 
reason why was already given by IStruvg (Il930l) : no Pi line is 
available in the "ordinary" range of the observed spectra of stars 
of spectral type F, G, or K. Some P ii and P iii lines are observ- 
able in the spectra of B-type stars. In the UV spectrum, some 
lines of Piv, Piii, and Pii can be detected in hot stars. To study 
the variation in the phosphorous abundance over time we need to 
use long-lived stars of type F, G, or K. In the era of the high res- 
olution spectrograph CRIRES- VLT, some lines of Pi, the lines 
of Mult. 1, at 1050-1082 nm, can be observed. For the first time. 
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it has thus become possible to trace the Galactic evolution of 
phosphorus. 

One single stable isotope of phosphorus, -^'P, is believed to 
be formed via neutron capture, as for the parent nuclei ^^Si and 
""'Si, probably in the carbon and neon burning shells during the 
late stages of the evolution of massive stars. The ""P produced in 
this way is the n released by the exp l osion of these massive stars 
as type II SNe. lWooslev & Weaverl(ll995h expect no significant 
production of phosphorus during the explosive phases. 

The handful of obs ervations available of P in diff'erent stars 
was summarised in Ca ff'au et al.l (12007). Since then. lHubrig et al.l 
(i2009i) have derived the P abundance in horizontal branch (HB) 
stars in globular clusters NGC 6397 and NGC 6752, where it 
appears to be strongly enhanced, by more than 2 dex and almost 
3 dex, respectively. Chemical anomalies in HB stars are usually 
attributed to d ifiiision effects, but the current mod els bv Michaud 
et al. (i2008i) do not seem to account for more than about 0.7 dex 
of the P enhancement in the photosphere of HB stars of similar 
temperature. 

Melendez et al. (|2009|) measured P in ten solar twins, al- 
though they do not state which lines they used. The lines of Mult. 
1 were not available in their observations and we speculate that 
they used the strongest line of Mult. 2 at 979.6 nm. 

A res ult concerning P evolu tion that we find intriguing, 
is that of ISbordone et al.l (|200^, who found that the sulphur 
abundance in the globular cluster 47 Tuc displays a signifi- 
cant correlation with the Na abundance. The authors point out 
that S could be produced through proton capture reactions on P 
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( P (p, y) S). However, P in the Sun is roughly 1.7 dex less Table 1. Infrared phosphorus lines analysed in this work, 
abundant than S. If the same abundance ratio held elsewhere, it 
would be impossible to produce any significant amount of S at 
the expense of P. The P abundance has never been measured in 
stars of 47 Tuc. 

In this paper, we use high-resolution infrared (IR) spectra, 
obtained with CRIRES at VLT, to analyse the behaviour of the 
P abundance as a function of me tallicity in the Galact ic disc. We 
selected a sample of stars from ICaffau et al.l ( l2005h for which 
the S abundance is known, to assess whether the proton capture 
on P is a viable production channel for S. The sample spans a 
metallicity range of 1.2 dex, -1.0 < [Fe/H] < +0.3. Our future 
plan is to extend the sample of metal-poor stars, and to include 
metallicities lower than [Fe/HQ =-1. 



X 


Transition 


J^low 


log 9f 


[nm] 




[eV] 




1051.1584 


4s'*Pi/2-4p*D0/2 


6.94 


-0.13 


1052.9522 


4s^P3/2^P^Dg/2 


6.95 


0.24 


1058.1569 


4s'*P5/2^P'*D0^2 


6.99 


0.45 


1059.6900 


4s^Pl/2^P*D?/2 


6.94 


-0.21 



2. Observed spectra 

The strongest of the Pi lines (1058. 1569 nm) has an equivalent 
width of about 2.2 pm in the spectrum of the solar photosphere. 
In a solar temperature star, all the lines of Mult. 1 are hardly 
visible when the metallicity is below [Fe/H] =-0.5. They are of 
high excitation energy (see Table[TJ and sensitive to temperature, 
hence become weaker in cooler stars. We wished to derive P 
abundance across as wide a range in metallicity as possible, to 
study the Galactic ev olution of phosphorus. From the sample of 
ICaffau et al.l ( l2005h . that is stars for which the abundance of 
sulphur was known, we selected 22 main-sequence or turn-off, 
F to early G-type stars. We gave preference to warmer stars, and 
tried to have a significant number of stars with detected planets. 
All the stars of the sample happen to be bright, in the range 3.3- 
7 mag in J. At the VLT-Antu 8 m telescope, these stars can be 
observed in twilight and under very poor weather conditions. 

All our spectra were observed by nodding along the slit. The 
slit width was 0'.'2, providing a spectral resolution R= 100 000; 
the adaptive optics correction was computed on axis, using the 
target star itself. Twenty of our stars have been observed in ser- 
vice mode in the program 386.D-0130 (PI, E. CafFau). We used 
the setting centred at 1059.5 nm in order 54, to ensure that the 
four detectable lines of P i were positioned on detectors 2 and 

3, which are the most efficient. The two reddest lines of Mult. 1 
(1068.1 and 1081.3nm), that are the weakest lines of the multi- 
plet in the solar spectrum, fall outside the range of the CRIRES 
detectors. Our analysis was performed with the reduced spectra 
provided by ESO. 



3. Model atmospheres and atomic data 

We analysed four of the IR P i line s of Mult. 1 . The atom ic data 
we used, listed in Table[Tl are from lBerzinsh et all (Il997h . 

For each star, we computed a ID-LTE model atmosphere 
with the code ATLAS 12 i n its Linux version (Kurucz 2005; 
ISbordone et al.l 120041 : ISbo rdone 2005), with the parameters 
(Teff/logg'/[Fe/H]) from ICaffau et al., (,2005.) and references 
therein, used to derive the sulphur abundance. The only excep- 
tion is HP 1461 fo r which we choose the atmospheric param- 
eters of ISousa et al. (2008) and the S abundance of Gonzalez 
Hernandez et al. pOlO) (see Sect.5 for details). This choice was 
driven by the need to make a consistent comparison between 
phosphorus and sulphur. The ATLAS 12 models were computed 
with a mixing length parameter of (Tmlt = 0.5 and the overshoot- 
ing option switched off. 



4. Analysis 

We measured the equivalent width (EW) of the P i lines, by us- 
ing the ira43task splot. The lines are weak: the strongest line 
has log(EW//i) - -5.36, but the majority of the lines have 
log(EW//l) < -5.5, and a Gaussian profile is a good approx- 
imation for the line profile fitting. The Pi line at 1051 nm is, 
for several stars, contaminated by telluric absorption (see upper- 
left panel of Fig.lTJ. When possible, we measured the EW of the 
P I line taking into account the presence of the telluric contam- 
ination, otherwise we rejected the line. Two of the lines (1052 
and 1058 nm) are blended with a Nil and Sii line, respectively. 
According to the strength of the lines and the stellar V sin /, we 
fitted the P i line alone, or took into account the close-by line at 
the same time, usi ng the deblending option of splot . We used 
the code WIDTH dKuniczll 19931 l2005l : ICasteUi! '2005') to derive 
the P abundance from the EW. The results are listed in Table|2] 
where the P abundance is given in the last but one column to- 
gether with the line-to-line scatter. Figure[T] shows the four P i 
Hnes for the star HD 13555. 

We looked at the effect that an uncertainty in temperature 
and gravity has on the P abundance determination. The compu- 
tations refer to the model 6112K/4.34/-H0.16. As a typical value 
for the uncertainty in the temperature and the gravity, we took 
ATeft = 100 K and Alog g = 0.2 dex, respectively. Looking at the 
literature, we found that Teff and log g in these stars usually do 
not differ between different authors by more than the uncertain- 
ties we assume here. The change in temperature ATeff = ±100 
induces A [P/H] =^gg^dex. For the gravity, Alogg - ±0.2 dex 
induces A [P/H] =^g jjj dex. The choice of the mixing length pa- 
rameter has very little influence on the determination of the P 
abundance. A change in the mixing length parameter, from 0.5 
to 1.5, increases the P abundance of +0.0045 dex, which is, for 
the quality of the data we analyse here, absolutely negligible. 
We assumed the systematic uncertainty in the P abundance to be 
0.1 dex, by linearly adding the contribution due to a change in 
effective temperature and in gravity. We neglect the contribution 
due to a change in the mixing length parameter. The statistical 
uncertainty due to the EW measurement is negligible with re- 
spect to this systematic error 

To our knowledge, no study of the deviation from local ther- 
mal equilibrium for phosphorus has been published, and no P 
model-atom, necessary for this c omputation, seems to exist in 
the literature. ICaffau et alJ (l2007h studied the effects that granu- 
lation has on the P abundance determination in the solar photo- 
sphere. These effects are very small (0.02, 0.03, 0.04, and 0.02 
for the four Pi lines we analyse here) compared to our error 
estimate, and we decided not to apply any three-dimensional 



[X/H]= log(Nx/NH) - log(Nx/NH)o and [X/Y]=[X/H]-[Y/H] 



' Image Reduction and Analysis Facility, written and supported 
by the IRAF programming group at the National Optical Astronomy 
Observatories (NOAO) in Tucson, Arizona, http://iraf.noao.edu/ 
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Table 2. Stellar parameters and abundances of our program stars. 



Target 



Teff 
K 



Logg [Fe/H] [S/H] 



Ref S/N 

#2/#3 



EW [pm] 
1051.1 



EW [pm] 
1052.9 



EW [pm] 
1058.1 



EW [pm] 
1059.6 



[P/H] Notes 



HD 1461 


5765 


4.38 


+0.19 


-0.05 


GIO 


180/150 




2.00 


2.70 


0.95 


+0.14 ±0.01 


P 


HD 13555 


6470 


3.90 


-0.27 


-0.25 


T02 


400/300 


0.80 


1.80 


2.90 


0.84 


-0.28 ± 0.03 




HD 25704 


5792 


4.20 


-0.91 


-0.71 


C05 


300/150 


0.30 


0.65 


1.20 


0.30 


-0.55 ± 0.05 




HD 33256 


6440 


3.99 


-0.37 


-0.30 


T02 


350/250 




1.90 


3.20 


1.00 


-0.20 ± 0.04 




HD 69897 


6227 


4.20 


-0.50 


-0.34 


C02 


120/50 


0.83 


2.00 


3.00 




-0.17 ±0.06 




HD 74156 


6112 


4.34 


+0.16 


-0.13 


E04 


250/90 


1.15 


2.40 


3.60 


1.00 


+0.05 ± 0.04 


P 


HD 75289 


6143 


4.42 


+0.28 


-0.03 


E04 


180/80 


1.50 


2.70 


3.60 




+0.16 + 0.01 


P 


HD 84117 


6167 


4.35 


-0.03 


-0.10 


E04 


150/70 




2.40 


4.00 




+0.09 ± 0.07 




HD 91324 


6123 


3.95 


-0.60 


-0.49 


C02 


300/140 




1.60 


2.40 




-0.30 + 0.02 




HD 94388 


6379 


3.96 


+0.07 


+0.12 


C02 


200/120 


1.40 


3.40 


4.40 




+0.03 ± 0.07 


V 


HD 120136 


6339 


4.19 


+0.23 


+0.05 


E04 


250/170 


1.90 


4.60 




1.90 


+0.26 ± 0.08 


P,V 


HD 139211 


6231 


4.12 


-0.26 


-0.16 


C02 


300/240 


1.40 


3.20 


3.50 


1.40 


+0.03 ± 0.05 




HD 207129 


5910 


4.42 


0.00 


-0.05 


E04 


400/250 


0.75 


1.60 


2.10 


0.70 


-0.09 ± 0.02 


T 


HD 207978 


6400 


4.03 


-0.63 


-0.53 


T02 


300/170 


0.69 


1.40 


1.60 




-0.47 ± 0.05 




HD 209458 


6117 


4.48 


+0.02 


-0.20 


E04 


200/100 




2.00 


3.00 




0.00 ± 0.03 


P 


HD 213240 


5984 


4.25 


+0.17 


-0.10 


E04 


300/180 


1.40 


2.75 


3.60 


1.20 


+0.15 + 0.01 


P 


HD 215648 


6158 


3.96 


-0.24 


-0.19 


C02 


220/100 


0.96 


2.10 


2.80 




-0.19 + 0.02 




HD 216385 


6300 


3.91 


-0.27 


-0.18 


T02 


400/180 


1.10 


2.40 


3.00 


0.90 


-0.18 + 0.02 




HD 216435 


5938 


4.12 


+0.24 


+0.10 


E04 


250/130 


1.80 


3.30 


4.10 


1.50 


+0.26 ± 0.02 


P 


HD 222368 


6178 


4.08 


-0.13 


-0.07 


C02 


280/130 


1.30 


2.40 


3.50 


1.20 


-0.03 ± 0.02 


V 




lef" refers to the reference for the stellar parameters and [S/H], and the 
lidai et al.U2002h: C05:ICaffau et alJ (12005b: C02:IChen et al.U2002h: E07 


corresponds 










The column " 


to, G10:lGonzalez 


Hernandez et alj (120101): 


T02:lTakada-F 


:lEcuvillonetal.l(l2004^. 
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Fig. 1. The four P i lines are here shown in the case of HD 13555. 
The observed spectrum (solid black) is compared to the synthetic 
profile (solid red) with the P abundance derived from the EW 
measurement. Also shown are the synthetic profiles obtained by 
changing the P abundance by ±0.1 dex (solid green). The tel- 
luric absorption affecting the 1051pm line is clearly visible in 
the upper-left panel. 



(hereafter 3D) correction. The stars of our sample span a small 
rang e in effective temper ature and gravity and our 3D-model 
grid dLudwig et al.ll20()9l) is very coarse. An application of the 



3D corrections interpolated from the small number of avail- 
able models would result in tiny shifts in the abundances, much 
smaller than the associated error of ±0.1 dex.. 



5. Discussion 

In Fig.|2| [P/Fe] is plotted versus the metallicity, [Fe/H]. The 
[P/Fe] ratio has a behaviour similar to that of an c-element. We 
can clearly see [P/Fe] increases as [Fe/H] decreases at metallic- 
ities as low as solar metallicity, [P/Fe] is close to zero for solar 
metallicity stars, and increases as metallicity decreases. 

For the light odd-Z elements it is generall y assumed that 
the m ain production channel is neutron capture (iBurbidge et al.l 
Il957h . If this is the case, one can expect that the production of 
P is sensitive to the neutron exces^ 77 . Since the neutron excess 
decreases with decreasing metallicity one can naively expect that 
the ratios of P to elements that are insensitive to 77, such as S or 
Fe, should decrease with decreasing metallicity. The increase in 
[P/Fe] with decreasing metallicity contradicts this simple expec- 
tation. It should however be noted that in the metallicity range 
covered by our study neither Na nor Al, the nearest odd-Z ele- 
ments, disp lay a decrease in the i r ratios to iron with decreasing 
metallicity dGehren et alj|2004l : lshret al.ll2004l : rGehren et alJ 
2006), but instead remain constant. 

Owing to the small number of metal-poor stars, it is unclear 
whether [P/Fe] reaches a plateau at about [Fe/H] a; -0.5, or con- 
tinues increasing. More data at metallicity around -1.0 would be 
very useful to clarify this interesting question. The observed data 
in Fi g.Elare compared to the model of P evolution in the Galaxy 
by Koba vashi et al.l ( 1200 6') (solid line). The model seems to qual- 
itatively reproduce the observed behaviour, although quantita- 



^ The neutron excess is defined by lAmetll ( Il97lh as: 77 = («„ - 
"p)/{n„ + «p), where n„ represents the total number of neutrons per unit 
mass and n^ the number of protons per unit mass. 
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Fig. 2. [P/Fe] as a function of the metallicity, [Fe/H]. The di- 
mension of the symbols reflects the average line-to-line scatter 
of 0.045 dex. The error bars are the sum under quadrature of 
the uncertainties of P (the linear sum of line-to-line scatter and 
the systematic uncertainty) and of Fe, the latter assumed to be 
0.1 dex for all stars. Stars known to have planets (red squares) 
have different symbols from stars without known planet (black 
hexag ons). A model of Galactic evolution of P ( Kobavashi et al. 
|2006|) is added (solid line) for comparison to the observed data. 



tively the increase in [P/Fe] is stronger than predicted and al- 
ready starts at solar metallicity. 

In our sample, seven stars have detected planets and are plot- 
ted as squares in Fig.|2l We see no difference in the behaviour 
in the [P/Fe] trend between the sample with and without plan- 
ets. One should always bear in mind that for the stars "without 
planets" it cannot be excluded that they harbour undetected plan- 
ets. We note, however, that all of the stars "without planets" in 
our sample, except for HD 215648 and HD 207978 have several 
high precision radial velocity points measured with the HARPS 
spectrograph at the ESO 3.6 m telescope. Given that, based on 
these radial velocity measurements, no planet has been reported, 
short-period massive planets ("hot Jupiters") should be excluded 
for these stars. 

In Fig. [3] [S/Fe] versus (vs.) [Fe/H] is shown for the same 
sample of stars. There is a tiny difference between stars with 
and without planets. The difference is reflected in the [P/S] vs. 
[Fe/H] plot in Fig.H] Phosphorus and sulphur behave similarly 
when effective temperature and gravity of the model atmosphere 
are changed. We computed the systematic uncertainties in S with 
the models we used for the P analysis, and assume that the sys- 
tematic error in the [P/S] ratio is given by the difference between 
the (logarithmic) systematic errors in P and S. The stars with 
planets are the most metal-rich ones in the sample, and have a 
lower sulphur abundances than those without planets. The trend 
of [P/S] with the metallicity is flat, and has a large scatter (see 
Fig.©. The average value of this "plateau" is [P/S]=0.10 ± 0.10. 
As noted above when discussing the [P/Fe] ratio, this "plateau" 
in [P/S] is surprising. In the same metallicity range, the ratios of 
either Na o r Al to the a eleme nt Mg decrease with decreasing 
metalUcity JGehren et alj|2006h . 

For the seven stars with known planets, we find ([P/S]) - 
0.15 ± 0.10, while for the stars without known planets ([P/S]) - 
0.07 ±0.10. This finding is not statistically significant, because 
the two values agree to within one cr. We also note that there 



Fig. 3. [S/Fe] as a function of the metallicity, [Fe/H]. The error 
bars indicate the uncertainties of [Fe/S] as given in the related 
papers. The symbols have the same meaning as in Fig.|2] 



is hardly any overlap in metallicity between the stars with and 
without planets. Nevertheless, this difference is puzzling and 
caused by the evolution in the abundance of sulphur, not phos- 
phorous. 

The two suspected binary stars (HD 33256, HD 69897) 
do not behave very differently from the complete sample. 
We can therefore conclude that the companion is either of 
the same spectral type or much fainter In addition, the 
three stars classified as variable in the SimbacQ astronomical 
database (HD 94388, HD 120136, HD 222368) and the pre main- 
sequence star (HD 207 129) do not display any difference from 
the other stars in the sample. For the three variable stars, the pho- 
to metric variab ility of HD 94388 and HD 222368 was reported 
by iPetili ( Il990l) . Later, however, these two stars were used as 
standards by several authors (e.g. HD 94388 for delta Scuti stars 
by Hintz et al. ( 1998), and HD 222368 for long-term variables 
by Manfroid et al. ( 1995). The shght vai-iabihty of HD 120136 
is studied by Hall et al. (2009), who confirm the long-term vari- 
ability detected by Lockwood et al. (2007). 

For HD 1461, we r eport in T able[2] the measurement of the 
sulphur abundance by Gonzalez Hernandez et al j ( 2010f) . who 
used the atmospheric pa rameters of lSousa et alJ (l2008h . 

Rivera et al. (1201 Ol) found a super-Earth on a short-period 
orbit (5.77 days) around HD1461. This star, and its planet, 
has been the subject of an intense debate between Gonzalez 
Hernandez et al. (2010) and .Ramirez et al. (20 10.). The de bate 
revolves around the claim made by|Melendez et alj ( l2009l) that 
the Sun has a peculiar chemical composition with respect to the 
average of solar-twins and solar analogs, being overabundant 
in volatile elements and under-abundant in refractory elements. 
iMelendez et alJ (l2009h proceed to identify the reason for this 
"peculiarity" is the presen ce of rocky planets, su ch as the Earth. 
This claim is supported bv'Ramir ez et alj ( 1201 Ol) . but rejected by 
Gonzalez Hernandez et al. (2010i) who, in an independent anal- 
ysis, do not confirm the results of IMelendez et al] (l2009l) and 
point out that HD 1461, hosting a super-Earth, has an abundance 
pattern that differs from what is cl aimed to be a signature of 
the presence of terrestrial planets by Melendez et alj 0200 9) and 
Ramirez et al. (201 0). This conclusion is contested by Ramirez et 



http://simbad.u-strasbg.fr/simbad/ 
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Fig. 4. [P/S] as a function of the metallicity, [Fe/H] . Symbols are 
as in Fig. 12] The error bars are computed as the linear sum of 
th(! systematic error of [P/S] (see text) and V2 times the P line- 
to- line sc; itter (assuming that random errors of P and S are of 
siiiilar si2e). 



al. (1201 Ol) . who claim, when restricting the analysis to refractory 
elements, that HD 1461 conforms to their claim. In this debate, 
the abundance of phosphorous provides data for a ne w element. 
In the sample of solar twins of iMelendez et al.l (l2009h . the mean 
[P/Fe] is -0.051 with a cr of 0.063, in excellent agreement with 
what we found in HD 1461 (-0.050), despite being rocky planet. 

6. Conclusions. 

We have determined the phosphorus abundance in twenty stars 
of the Galactic disk that we observed with CRIRES at VLT. The 
metallicity of the sample extends over about 1 .2 dex, providing a 
first insight into the evolution of phosphorus in the Galaxy. The 
only published model of the chemical evolution of the Galaxy 
that pr ovides informatio n about the evolution of phosphorus is 
that of iKobavashi et all (I200&) . The model predicts correctly, 
from a qualitative point of view, the rise of the P/Fe ratio with 
decreasing metallicity. The observed rise is, however, larger by 
almost a factor of two than the model's prediction, and lacks the 
"plateau"' at metallicities [Fe/H] > -0.5. Our analysis indicates 
that phosphorus behaves in a similar way to an a-element, such 
as sulphur. The model of Galactic chemical evolution described 
by Gibs on et al.l ([2005) predicts a ~ 0.2 dex enhancement of the 
[P/Fe] ratio over two dex in metallicity, from solar to -2.0. The 
level of enhancement at low metallicity is compatible with what 
we observe, but no mention is made of the decrease in this ratio 
near solar metallicity. 

We have found that the [P/S] is roughly constant at a level of 
about 0. 1 dex. However, the star-to-star scatter makes this value 
compatible with zero. This constant value implies that P and S 
are produced in the same amounts over the metallicity range con- 
sidered by our investigation. The picture is different from what 
is found for Na or Al: both [Na/Mg] and [Al/Mg] decrease sig- 
nificantly with decreasing metallicity. A possible interpretation 
is that the P production is insensitive to the neutron excess 77. If 
this is the case, then P should be produced by nuclear reactions 
other than neutron captures. Possibilities that should be explored 
are proton captures on ^"Si and a captures on ^^Al. A corollary 
of this analysis is that, since in the Sun sulphur is about a fac- 



tor of 30 more abundant than phosphorus, this ratio remains the 
same at all metallicities (in the range that we have examined). 
We therefore conclude that S production through proton cap- 
tures on P is not a viable mechanism. Although the possibility 
of a large enhancement of P abundance in globular clusters can- 
not be ruled out, it seems unlikely. 

It is interesting to compare our results to phosphorus abun- 
dances derived from measurements in damped Ly-a (DLA) sys- 
tems at high redshift. These are gas-rich galaxies observed as 
absorption lines in the direction of QSOs. iMolaro et al.l (l2001h 
measured the P abundance in the DLA at z = 3.3901 towards 
QSO 0000-2620. This system has a very low metallicity and 
measurements of [P/H]=-2.31, [P/Fe]=-0.27, and [P/S]=-0.40. 
The ratios of P to iron and sulphur are clearly different from 
what those seen in Galactic disk stars. This may be due to the 
low metallicity of the system, suggesting that perhaps also in 
our Galaxy, the [P/Fe] ratio starts at low values, increases un- 
til reaching a maximum, and thereafter decreases towards so- 
lar metallicity. This be haviour is predicted by the model of 
IKobavashi et al.l (l2006t) . although not shown in Fig .2, which 
only extends down to [Fe/H] =-1.0. Another possibility is that 
this DLA galaxy has undergone a different chemical evolution, 
since it has [S/Fe] =-H0.13, while at th ese metallicities [S/Fe]~ 
-I-0.4 in the Galaxy dCaffau et al.l l2005h. 

Fenner et al. (|2004|) determined an upper limit of [P/S]< 
-hO.01 in the DLA at z = 2.626 towards FJ081240.6-H320808. 
This system is more metal-rich ([O/H]=-0.44), thus probably 
of metallicity comparable to some of the stars in our sample. 
However, the [P/S] ratio is slightly lower than what we observe. 
Although this DLA system may resemble the Milky Way disk, it 
has probably experienced a different chemical evolution. 

Our investigation is the first step towards an understanding 
of the evolution of P in the Galaxy. It will be important to secure 
more data at metallicity [Fe/H] =-1.0 and lower, to understand 
whether the [P/Fe] ratio flattens out, continues to rise, or finally 
decreases towards low metallicity. The good performance of the 
CRIRES spectrograph suggests that it will be possible to observe 
considerably fainter targets with adequate exposure times. A di- 
rect measurement of P in 47 Tuc or other globular clusters would 
allow us to probe the feasibility of S production through proton 
capture in these objects. 
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